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ON AN ELLIPTIC EQUATION ARISING EROM PHOTO-ACOUSTIC 
IMAGING IN INHOMOGENEOUS MEDIA 

HABIB AMMARI, HONGJIE DONG, HYEONBAE KANG, AND SEICK KIM 


Abstract. We study an elliptic equation with measurable coefficients arising from 
photo-acoustic imaging in inhomogeneous media. We establish Holder continuity 
of weak solutions and obtain pointwise bounds for Green's functions subject to 
Dirichlet or Neumann condition. 


1. Introduction and main results 
We consider an operator L of the form 

Lu = - drv(y(x)Vu) + iku 

acting on a complex-valued function u defined on a domain Q c R" (n > 3). Here, 
i= V=I, A: is a nonzero constant, and y(x) is a real n X n matrix valued function 
satisfying 

v|£| 2 <R eFy(x)Z, e C", VieO; ||ylk~ (n) < v~\ (1.1) 

for some v e (0,1]. The adjoint operator L* is given by 

L*u = - div(y(x) T Vu) - iku. 

Note that the condition lll.lt is invariant if we replace y by y T . Therefore, even if 
L is not necessarily a symmetric operator, L and L* are operators of the same type. 

The operator L arises from quantitative photo-acoustic imaging in inhomoge¬ 
neous media; see HUH- The purpose of quantitative photo-acoustic imaging is 
to reconstruct the optical absorption coefficient of an optical anomaly from the 
absorbed energy density. The absorbed energy density is obtained from boundary 
measurements of the pressure induced by the thermo-elastic effect; see (10). The 
reconstruction problem is challenging. One promising method is implemented in 
El in the framework of diametrically small optical anomalies. It uses asymptotic 
expansions of the energy density in terms of the characteristic size of the optical 
anomalies. Such expansions are based on estimates of the Neumann function of 
the operator L under the assumption that y is uniformly continuous in Q. This 
note aims at removing such regularity assumption and extending the results of |U 
to real biological tissues. 
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We say that the function u e H 1 (D) is a weak solution of Lu = f, f e L 2n/( - n+2 \Q), 


if it satisfies for any <p e Hj(Q) the identity 


| (V<p T yVu + ikucp) = | fcp. 
Jo Jo 


( 1 . 2 ) 


In this short article, we prove that any weak solution u of Lu = 0 satisfies the 
following De Giorgi-Moser-Nash type estimate in B r (xf) c Q: 

[u\a;B r/2 (xo) < Cr“ n/2 “ a ||u|| L2(Br(lo)) , 

where [u] a -Bj Xo -) is the standard Holder semi-norm of u over an open ball of radius 
r centered at Xq. We emphasize that the constant C in the above estimate is inde¬ 
pendent of k and depends only on n and v. Since u is complex valued and its real 
and imaginary parts are intertwined in the equation through the parameter k, this 
result is not immediate from the classical De Giorgi-Moser-Nash theory. We also 
establish similar estimates near the boundary for weak solutions satisfying zero 
Dirichlet or Neumann boundary condition. Then by utilizing these estimates, we 
construct Dirichlet and Neumann functions and derive the pointwise estimates for 
them; we show that they are bounded in modulus by C|x - i/| 2- ". By Dirichlet's 
function we mean a function G(x, y) satisfying 

LG(-, y) = by in Q, 

G(-, y) = 0 on <9Q, 

and similarly, by Neumann function, we mean a function N(x, y) satisfying 

LN(-, y) = by in Q, 
yVN(-, y) ■ n = 0 on dQ. 

The above definitions for Dirichlet and Neumann functions should be interpreted 
in a weak sense; see 0 Sec. 2.2] for details. 

Now, let us state our main results more precisely. 

Theorem 1.3. Let u e H 1 ( Q) be a weak sohition ofLu = f. 

(i) Let B r = B r (x o) c Q and f e U\B r )for some p > «/2. Then, we have 

[u\ a;Br/2 < Cr~ a (r - ” /2 ||i;|| L 2 (Br) + r 2 ^ ||/|b (Br) ), (1.4) 

where a e (0,1) and C > 0 are constants depending only on n, v, and p. 

(ii) Suppose u = 0 on 30, in the sense that u e Hj(Q). Let Q r = Q r (xo) = ^ n B r (xo) 
with xq e Q and f e LP{D. r )for some p > n/2. Then, we have 


IMlL“(Q r/2 ) <c(r " / 2 ||w||l2(d,) + r 1 " /p ll/lliP(Q r )) 


(1.5) 


for any 0 < r < diam Q, where C is a constant depending only on n, v, and p. 

(Hi) Suppose that Q is a bounded Lipschitz domain and du/dv = 0 on <90 in the sense that 
for any cp e H 1 ( Q) the identity (IT.2b holds. Then, there is a constant C depending 
only on n, v, p, and Q such that the estimate 111.5b is valid for any 0 < r < diam Q. 

Remark 1.6. It should be clear from the proof that the same conclusions of the 
preceding theorem hold true for any weak solution u e H 1 (LI) of L*u = f. 


ELLIPTIC EQUATION ARISING FROM PHOTO-ACOUSTIC IMAGING 


3 


Remark 1.7. In view of the proof of Theorem 11.31 below, by using the boundary 
De Giorgi-Moser-Nash theorem, it is easily seen that under the assumption of 
Theorem 11.31 i i) and a uniform exterior cone condition on Q, we have the following 
Holder estimate: 

[«W < Cr~ a (r~ n/2 \\u\\L 2 (Q,) + t 2 -” /p ll/lb’(Q,))' 

where a e (0,1) and C > 0 are constants depending only on n, v, p and the 
parameters in a uniform exterior cone condition. The same estimate also holds 
under the assumption of Theorem II ,3l iii). 

Theorem 1.8. (Dirichlet's function) There exists a unique Dirichlet function G(x, y) such 
that 

| G(x, y) | < C|x - y\ 2 ~ n , Vr,j/eQ, x + y, (1.9) 

where C is a constant depending only on n and v. In particidar, when Q = R", we have 
the fundamental solution with the pointwise estimate & 

Theorem 1.10. (Neumann function) Let Cl be a bounded Lipschitz domain. There exists 
a unique Neumann function N(x, y) such that 

\N(x,y)\ < Cmax(l, \k\~ 1/2 ) \x- y\ 2 ~ n , Vx, y eO, x + y, (1.11) 

where C is a constant depending only on n, v, and Cl. 

We note that a similar result regarding Neumann function N(x, y) was estab¬ 
lished in a under the stronger assumption that y is uniformly continuous in Q 
and Q is a C 1 -domain. Therefore, our result is quite an improvement of that in f4j. 
We emphasize, however, that in [S4j estimates of derivatives of the Neumann func¬ 
tions are also derived under additional smoothness assumptions on the coefficients 
and the domain. 

Finally, one may wish to consider the case when k is not a constant but a variable 
coefficients k(x) and include lower order terms in L. Even in this case, our main 
estimates are independent of the size of k(x) but on the size of its "oscillation". We 
shall illustrate how to generalize our results to this case in Section [Z4l 

2. Proof of main theorems 

2.1. Proof of Theorem 11.31 The idea of the proof is to add a new variable and 
rewrite the equation as a parabolic equation without zero-order term, so that k 
does not play a role in our estimates. Let v(x, t) = e ,kt u(x). Then v becomes a weak 
solution of the parabolic equation 

v t - div(y(x)Vv) = e lkt f(x) in Ox (- 00 , 00 ). 

Note that the real and imaginary parts of v are solutions of a real scalar equation 

w t - div(y(x)Vrp) = g(x, t) in Q x (- 00 , 00 ) 

with g(x, t) = cos (kt)f(x) or g(x, t) = sin(kt)f(x). Denote Q r = B r x (-r 2 , 0). Clearly, 
we have 

^ WfWu’iBr) 

and thus by the De Giorgi-Moser-Nash theorem applied to the real and imaginary 
parts of v, we have the parabolic Holder estimate 

[»W(Q, /2 ) ^ Cr-(r- ( - 2)/2 ||u|| L2 ( Qr)+ r 2 -^||gll Lr L;( B ,))^ 
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where C depends only on n, v, and p. 

Then for any x, y e B r / 2 (^ 0 ) with x + y, we have 

I u(x) - u(y)\ _ \v(x, t ) - v{y, t)\ 

\x - y\ a \x - y\ a 

< Cr~ a (r-< n+2 ) /2 \\v\y lQr) + r 2 -’ ,/p llglli. ;r(Qr) ) 

= Cr~ a (r~ n ^ 2 \\u\\L 2 (B r ) + r 2 ~ n/r \\f\\ip(B r )j ■ 

Therefore, we have the estimate m. This completes the proof for part (i). 

Next suppose that u e Then v(x, t ) = e lkt u(x) vanishes on <9Q x (- 00 , 00 ). 

Denote U r = O r X (—r 2 , 0). By the local maximum principle of De Giorgi and Moser 
(see, e.g., 0 Ch. Ill, § 7]) applied to both real and imaginary parts of v, we get 

IMIl“ ( u, / 2 ) < c(r-(" +2 )/ 2 ||u|| L2(Ur) + ^llgll^tu,)), (2.1) 

where C depends only on n, v, and p. Then as in part (i), we get (11.5b from the above 
inequality. This completes the proof for part (ii). The proof for part (iii) is similar. 
The assumption that Q is a Lipschitz domain makes the (multiplicative) Sobolev 
inequality available to us and we can still apply De Girogi or Moser's method as 
mentioned above to get the local boundedness estimate Il2.1b . The constant C at 
this time, however, additionally depends on the constant of the Sobolev inequality 
and diameter of Q, and thus depends on Q; see |j5) § 5.3.1]. ■ 


2.2. Proof of Theorem ll.8l First, we establish unique solvability of the problem 

Lu=f in Q, / e L 2n/( " +2) (Q) (2.2) 

in the space H k ( Q). For any R> 0, denote Qr = Q n Br(0) and consider the Hilbert 
space Hq(Qr) equipped with inner product 

(u, v) = I Vu ■ Vzi. 

Jo,r 


We define the sesquilinear form B : H k ( Or) X H 2 (D.r) —» C as 


B[u, v] 



yVu + iknv). 


It is easy to check that B satisfies the following: 
i) there exists C = C(n, v, k, R) such that 

\B[u,v]\ < C||m|||M|, Vu,ueHj(Q R ); 


ii) 


r'||n|[ 2 < ReB[u,u], Vh e Hq(Qr). 
Let F : Hq(Qr) —» C be a bounded anti-linear functional; i.e., 

F(v + Aw) = F(v) + AF(w), Vu, w e Hq(Qr), 


VAeC 


and 

l|f||* := sup|F(u)| < 00 . 
IWI<1 
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(2.4) 


(2.5) 


Then by the Lax-Milgram lemma (in the complex case), there exists a unique 
u E Hj( Qr) such that ||u|| < v _1 ||F||, and 

B[u,v] = F(v), VveHl(Q R ). 

In particular, for any / E £ 2 «/(»+ 2 )( q ), there is a unique u <K> E H*(Q R ) such that 

f (Vv T yVu (R) + iku (R) v) = f fv, Vv e Hj(Q R ) (2.3) 

CIr J Qr 

and «(*) satisfies the estimate 

I|Vw (R ^Hl 2 (q r ) < C(n)v 1 11/1 |i, 2 n/(n+ 2 )( 0 ). 

By taking i> = in d2.3b and using Cauchy's inequality and the Sobolev embed¬ 
ding theorem, we get 

Jw llu (R) llu 2 (o R ) ^ C(n,v)||/|| L 2 »/<n+ 2 )( Q ). 

Therefore, there exist a sequence (My] and a function u in Hj(Q) such that My —>■ u 
weakly in H 1 (D), and m satisfies 

I|Vm|| L 2 (q) < C(n)v 1 ||/|| L 2 1 ./( 11 + 2 )(Q), yf\k\ ||m|| L 2 (q) < C(n,v)||/|| L 2 n/(n+ 2 )(Q). 

Thus by taking the limit in d2.3b . we have 

I (Vv T yVu + ikuv) = f fv, VneC“(Q). 

Jn J a 

By d2.4b and a density argument, we get d2.5b for any v E JTj(Q). This establishes 
existence of a solution in the space Hj(Q). In the case when / = 0, the above 
identity implies that if m e Hj(Q) is a weak solution of Lu = 0, then u = 0. This 
establishes uniqueness too. 

Having unique solvability of the problem d2.2b and the part i) and ii) of The¬ 
orem [L3l one can construct Dirichlet's function and obtain the pointwise bound 
ill.9b by following the same method as in 0 and 0- The uniqueness of Dirichlet's 
function is a mere consequence of the fact that u = 0 is a unique weak solution of 
Lu = 0 in the function class of Hj(Q). ■ 

2.3. Proof of Theorem ll.101 The theorem follows from replicating the proof of fU 

Theorem 2.3] using part (i) and (iii) of Theorem ll.3l It is assumed that k > kg > 0 in 
0. However, by replacing the constant C = C(n, v, kf in the estimate 0 (2.29)] by 
max(l, \k\-^ 2 ) C(n, v), we dispense with introducing the lower bound /c () . ■ 

2.4. A generalization. In this subsection, we show how to extend our results to 
operators with lower-order terms when k is not a constant but a variable coefficients 
k(x). We consider an operator L of the form 

Lu = - div(y(x)V« + b(x)u) + b(x) ■ Vm + c(x)u + ik(x)u, (2.6) 

where b = (b 1 , ..., b”) T , b = (b 1 ,..., b n ) T , and k(x) are real valued functions defined 
onQclR". We assume that b, b, and c are bounded functions so that 

Pllco + ll&lloo + Hello. < M. 

When k(x) = /c e IR, by the same proof, the estimate dl.4b in llieorem l l .31 (1) remains 
valid for all r E (0, r 0 \, where r 0 > 0 is any fixed number, but the constant C now 
depends on v, n, p, r 0/ and M; but it is still independent of k. 
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Next we treat the case with general variable coefficient k(x). We assume that k(x) 
satisfies 


sup r 2 ( -r | k(x) - k r | I? dx\ < k 0 , 
re(0,r o ] \JB r I 


where 


k r = 



dx, 


for some r 0 > 0, q > n/2, and k„ > 0. This condition is satisfied, for instance, when 
k(x) is uniformly continuous or is a function of bounded mean oscillations (BMO). 
We rewrite the equation Ln = f into 

- div(y(x)V«) + b(x) ■ V» + c(x)u + ik r u = f - i(k(x) - k r )u. 

Choose pi e (n/2, q) and p 2 e (n/2, oo) such that 1/p! = l/q + l/p 2 . By the previous 
statement regarding an operator L with k(x) = k e ]R, we have 

[u\a;B rll < Cr~ a (r-” /2 || M || L2(Br) + /-n'nfWiHB,) + r 2 ~ n ^\\(k - kr)u\\m ( B r )) 

< Cr~ a (r- n/2 \\u\\ L 2 ( Br) + r 1 ~ nlv \\f\\u>(B r) + K 0 r~ n, P 2 \\u\\ m{Br )), 

where we used Holder's inequality in the second line. Estimate (11.4b then follows 
from the inequality above and a standard iteration argument (cf. [6} pp. 80-82]). 
Of course, now the constant C also depends on q and k 0 in addition to v, n, p, r 0 , 
and M, but is still independent of the size of k(x). 

Similar adjustments are needed for Theorem 11.31 (ii) and (iii). Theorem 11.81 
remains valid when Q is bounded and we have coercivity for the sequilinear form 
associated with d2.6b : i.e., under these assumptions, there exists a unique Dirichlet 
function G(x, y) such that 

I G(x, y )| < C\x - y| 2 “", 

where C is a constant depending on v, n, M, k„, and diam Q. We have the required 
coercivity, for example, when c(x) > c 0 for some c 0 large enough. Under an addi¬ 
tional assumption that k(x) is bounded away from zero. Theorem !! ,1 Ol also remains 
true if the estimate ill. lit is replaced by 

|N(x, y)| < Cmax (l, (inf \k\)~ 1/2 ^j \x - y\ 2 ~ n , 

where the constant C depends additionally on M and K 0 . 


Remark 2.7. As in the classical De Giorgi-Nash-Moser theory, it is possible to relax 
assumptions on b, b, and c so that they belong to certain Lebesgue or Morrey spaces 
other than the L°° space. We leave the details to the interested reader. 
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